dAM*  _ 

AIAA  97-0264 

Numerical  Investigation  of  Twin-Nozzle 
Rocket  Plume  Phenomenology 

Houshang  B.  Ebrahimi 
Sverdrup  Technology,  Inc.,  AEDC  Group 
Arnold  Engineering  Development  Center 
Arnold  Air  Force  Base,  TN  37389 

Jay  Levine 

Phillips  Laboratory/RKS 
Edwards  Air  Force  Base,  CA  93524 
and 

Alan  Kawasaki 

Sparta  Inc.,  Phillips  Laboratory/RKS 
Edwards  Air  Force  Base,  CA  93524 

CO 


35th  Aerospace  Sciences 
Meeting  and  Exhibit 

January  6-10,  1997  /  Reno,  NV 


For  permission  to  copy  or  republish,  contact  the  American  Institute  of  Aeronautics  and  Astronautics 
1801  Alexander  Bell  Drive,  Suite  500,  Reston,  VA  22091 


DISTRIBUTION  STATEMENT  A 

Approved  for  public  release; 
Distribution  Unlimited 


DTIC  QUALITY  mSPKOTED  0 


NUMERICAL  INVESTIGATION  OF  TWIN-NOZZLE  ROCKET  PLUME 

PHENOMENOLOGY* 

Houshang  B.  Ebrahimi,  **  Jay  Levine f  and  Alan  KawasakP 


Abstract 

The  Generalized  Implicit  Flow  Solver  (GIFS) 
computer  program,  developed  under  sponsorship 
of  the  Air  Force  Phillips  Laboratory,  Edwards  Air 
Force  Base,  CA,  has  been  modified  and  used  for 
three-dimensional  reacting  two-phase  flow  prob¬ 
lems.  The  intent  of  the  original  GIFS  development 
effort  was  to  provide  the  JANNAF  community  with 
a  standard  computational  methodology  to  simulate 
multiple  nozzle/plume  flow-field  phenomena  and 
other  three-dimensional  effects.  Recent  develop¬ 
ment  efforts  have  concentrated  on  improving  the 
run  time  and  robustness  of  the  algorithm. 

The  GIFS  computer  program  was  originally 
released  as  an  untested  research  version.  Since 
that  time,  several  corrections  and  enhancements 
have  been  made  to  the  model.  The  Van  Leer  Flux 
Splitting  option  has  been  successfully  implemented 
into  the  existing  GIFS  model  and  provides  a  more 
robust  solution  scheme.  A  Parabolized  Navier- 
Stokes  (PNS)  version  of  the  GIFS  algorithm  is  cur¬ 
rently  under  development  and  is  intended  to  sub¬ 
stantially  improve  the  run-time  requirements  for 
flow  fields  dominated  by  supersonic  flow  regimes. 
These  improvements  and  enhancements  will  foster 
the  application  of  the  GIFS  model  in  the  CFD  com¬ 
munity.  This  paper  reports  the  significant  results  of 
several  twin-nozzle/plume  applications  of  the  GIFS 
code.  Six  simulations  of  Titan  II  plume  flow  fields 
have  been  completed  to  assess  the  effects  of 
three-dimensionality,  turbulent  viscosity,  afterburn¬ 
ing,  near-field  shock  structure,  finite-rate  kinetic 
chemistry,  intranozzle  geometric  spacing,  and  ini¬ 
tial  nozzle  exit  plane  profile  effects  on  the  subse¬ 
quent  plume  exhaust  flow  field.  The  results  of 


these  calculations  indicate  that  the  viscous  stress 
model,  kinetic  chemistry,  and  nozzle  exit  profile  are 
significant  parameters  that  should  be  considered  in 
analyses  and  interpretation  of  the  calculations.  In 
addition,  the  intranozzle  geometric  spacing  influ¬ 
ences  the  Mach  reflection  location  which  can  sig¬ 
nificantly  affect  the  plume/plume  impingement 
shock  location,  inviscid  shock  structure,  and  shear 
layer  growth. 

Introduction 

In  order  to  support  testing  and  analysis  require¬ 
ments  of  the  plume  community,  a  need  exists  for  a 
fluid  dynamics  model  which  solves  the  fully  cou¬ 
pled  two-phase  Navier-Stokes  equations  in  multi¬ 
ple  dimensions.  Evaluation  of  solid-propellant 
rocket  motor  performance,  nozzle  erosion,  and 
solid-propellant  rocket  plume  radiative  transfer 
analyses  requires  a  computer  model  which  simu¬ 
lates  complex  three-dimensional,  chemically  react¬ 
ing  two-phase  flow  effects.1  Although  this  type  of 
full  Navier-Stokes  method  provides  an  accurate 
qualitative  description  of  the  basic  features  of  the 
propulsion-generated  flow  fields,  quantitative  simu¬ 
lations  for  predicting  fundamental  parameters  such 
as  base  pressure,  static  pressure,  temperature, 
and  chemical  composition  in  the  flow-field  domain 
have  not  been  validated.  In  the  past  few  years,  sig¬ 
nificant  progress  has  been  made  in  the  areas  of 
numerical  rocket  flow  simulations  and  computa¬ 
tional  resources  to  the  point  that  Navier-Stokes 
solutions  are  viable  analysis  tools. 

The  flow  fields  generated  by  rocket  propulsion 
systems  are  complex,  with  regions  of  strong  invis- 
cid/viscous  interactions,  free-stream  shear  layers, 
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nozzle  wall  and  missile  body  boundary  layers, 
external  and  internal  shocks,  separation  regions, 
and  plume/plume  impingement  and  associated 
flow  interactions  for  multiple  nozzle  designs,  all 
with  chemically  reacting  kinetics.2  Modeling  all  of 
these  phenomena  simultaneously  is  a  challenging 
task  and  is  pushing  the  state  of  the  art  (SOA)  for 
CFD  models.  To  account  for  all  phenomena  affect¬ 
ing  plume  flow  properties  and  the  resulting  radia¬ 
tive  transfer  implications,  computationally  efficient 
multidimensional  computer  models  are  required. 

In  the  recent  past,  there  has  been  increased 
emphasis  on  the  application  of  existing  CFD  mod¬ 
els,  both  in  the  commercial  arena  and  government- 
developed  computer  programs,  to  simulate  these 
complex  multidimensional  plume  phenomena. 
Almost  all  of  the  solutions  obtained  to  date  have 
used  the  perfect  gas  approximation,  as  fully  react¬ 
ing  flows  are  considerably  more  difficult  to  solve. 
Including  the  effects  of  chemistry  in  the  solution 
produces  a  stiff  set  of  equations  which  are  numeri¬ 
cally  difficult  to  solve  using  conventional  algo¬ 
rithms.  In  addition,  the  grid  resolution  requirements 
become  more  severe  when  reacting  chemistry  is 
included.  When  chemical  reactions  were  first  intro¬ 
duced  into  the  GIFS  Titan  simulation,  the  original 
version  of  the  GIFS  code  was  unstable  and  did  not 
produce  a  converged  solution.  This  was  a  difficult 
issue  to  resolve  and  required  a  special  enhance¬ 
ment  to  the  GIFS  algorithm. 

Conventional  rocket  exhaust  flow  computer 
models3,4  employ  Euler  solutions  for  the  plume 
core  flow  and  superimpose  a  turbulent  parabolic 
mixing  approach  for  the  air/plume  entrainment 
(shear)  layer.  These  methods  are  generally  not 
adequate  for  situations  when  the  flow  is  not  fully 
dominated  by  either  inviscid  core  expansion  or  the 
plume  afterburning  phenomena.  Three-dimen¬ 
sional  features  which  result  from  multiple  nozzle 
systems  and  vehicle/base  interactions  cannot  be 
sufficiently  treated  through  the  use  of  these  mod¬ 
els.  Inaccurate  accounting  of  the  upstream  influ¬ 
ences  in  the  plume  shear  layer  is  one  of  the  primary 
reasons  which  prevents  the  approximate  models 
from  accurate  simulation  of  the  overall  flow  field.5 

The  GIFS  numerical  algorithm  provides  a  solu¬ 
tion  of  the  two-  and  three-dimensional  Reynolds- 


averaged  Navier-Stokes  (NS)  equations  using  the 
MacCormack  implicit  finite-volume  algorithm  with 
Gauss-Seidel  line  relaxation.6  Several  2-D  and  3-D 
plume  flow-field  calculations  have  been  completed 
for  the  plume  nearfield  region  using  the  original 
GIFS  model.6  The  GIFS  model  includes  a  frozen 
and  a  generalized  finite-rate  kinetic  chemistry 
model,  a  Lagrangian  particle  model  for  treating 
solid  or  liquid  particulates,  and  a  two-equation  tur¬ 
bulence  model,  as  well  as  a  laminar  model.  These 
complex  phenomena  are  required  to  accurately 
simulate  the  physics  expected  to  contribute  to  the 
subsequent  plume  signature.  The  original  objective 
of  the  GIFS  development  was  to  provide  a  flow- 
field  model  to  be  used  in  conjunction  with  radiative 
transfer  models  to  predict  base  heating  character¬ 
istics  resulting  from  three-dimensional,  multiple 
nozzle  propulsion  configurations.  The  Van  Leer 
Flux  Splitting  option7  has  been  successfully  imple¬ 
mented  into  the  existing  GIFS  model  and  provides 
a  more  robust  solution  scheme.8  A  Parabolized 
Navier-Stokes  (PNS)  version  of  the  GIFS  algorithm 
is  currently  under  development  and  is  intended  to 
substantially  improve  the  run-time  requirements  for 
flow  fields  dominated  by  supersonic  flow  regimes. 
These  improvements  and  enhancements  will  foster 
the  application  of  the  GIFS  model  in  the  CFD 
community. 

The  plume  structure  of  a  multinozzle  hypersonic 
vehicle  flying  at  high  altitude  is  largely  dominated 
by  the  processes  taking  place  in  the  plumes’  inter¬ 
action  region.  Prior  to  the  release  of  the  GIFS 
model,  multiple  nozzle/plume  flow  fields  were  com¬ 
monly  treated  by  assuming  a  single  equivalent  noz¬ 
zle  configuration  having  equal  mass,  energy  and 
momentum  of  the  multiple  nozzle  geometry.  Fur¬ 
ther,  a  uniform  nozzle  exit  flow  was  used  as  the 
starting  conditions  for  the  plume  calculation.  The 
simplified  model  assumes  that  the  details  of  the  3- 
D  flow  structure  in  the  near  field  flow  are  unimpor¬ 
tant  and  that  the  flow  processes  which  affect  the 
plume  shear  layer  initialization  (such  as  base  sepa¬ 
ration  and  recirculation)  will  be  dominated  by  the 
overall  ambient  flow  entrainment  effects.  The  level 
of  agreement  between  computations  using  the  sin¬ 
gle  equivalent  nozzle  methodology  and  radiometric 
data  from  multiple  nozzle  propulsion  systems  has 
not  been  acceptable,  especially  for  spatial  analysis 
requirements.  The  source  of  the  disagreement  is 
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due  at  least  in  part,  to  an  incorrect  physical  model 
of  the  phenomena  dominating  the  observations  and 
can  be  further  identified  in  the  model  as  incorrect 
turbulence,  chemistry,  and  an  inadequate  repre¬ 
sentation  of  the  geometry. 

The  motivation  for  this  study  is  to  demonstrate 
the  significance  of  three-dimensional  effects  as 
applied  to  multiple  nozzle  rocket  missile  plumes,  to 
determine  how  these  phenomena  may  be  better 
simplified  in  order  to  promote  improvements  to 
engineering  approaches,  and  to  explore  methods 
to  reduce  the  overall  CPU  resource  requirements 
for  simulating  three-dimensional  solutions.  The 
study  presented  in  this  paper  represents  the  initial 
efforts  in  simulating  the  rocket  exhaust  flow  field 
from  a  hypersonic  vehicle  in  flight. 

Computations 

The  computational  effort  consisted  of  six,  three- 
dimensional  twin-nozzle  calculations  for  the  Titan  II 
vehicle  at  flight  conditions.  In  all  cases,  only  the 
plume  flow  field  was  computed,  i.e.,  base  effects 
were  not  considered.  (Complete  solutions  including 
the  flow  over  the  missile  body,  the  base  region, 
and  the  exhaust  plume  domain  will  be  presented  in 
a  future  paper.)  Also  in  the  first  five  cases,  bound¬ 
ary-layer  effects  from  the  rocket  nozzle  were  not 
considered.  The  starting  boundary  conditions  were 
assumed  to  be  uniform.  In  the  last  case,  the  plume 
was  computed  utilizing  two-dimensional  radial  pro¬ 
file  starting  conditions  at  the  exit  plane  location  of 
the  nozzle.  These  were  provided  using  the  Two- 
Dimensional  Kinetics  Computer  Program,  TDK.9 
For  all  cases,  the  plume  flow  field  was  simulated  at 
an  altitude  of  approximately  50  km  at  a  free-stream 
Mach  number  of  5.7.  The  resulting  nozzle  exit-to- 
free-stream  velocity  ratio  was  approximately  5,  and 


the  chamber  pressure  to  freestream  pressure  ratio 
was  approximately  50,000:1.  The  nozzle  operated 
in  all  cases  in  a  significantly  underexpanded  mode. 
The  parametric  flow-field  conditions  for  all  six 
cases  are  summarized  in  Table  1. 

The  Titan  II  Space-Launched  Vehicle  (SLV) 
propulsion  system  is  propelled  by  two  engines 
located  on  either  side  of  a  plane  of  symmetry  pass¬ 
ing  through  the  vehicle  centerline.  To  limit  the  num¬ 
ber  of  grid  points  and  CPU  run  time,  quarter  sym¬ 
metry  assumptions  were  made  for  the  computa¬ 
tional  domain  plane.  Zero  angle  of  attack  and  col- 
linear  plume  exhaust  were  two  additional  con¬ 
straints  imposed  for  these  simulations.  A  total  of 
600,000  grid  points  were  utilized  in  the  computa¬ 
tional  domain  as  shown  in  Fig.  la  with  the  orienta¬ 
tion  indicated  by  the  axes.  A  schematic  of  the  3-D 
plume  configuration  is  shown  in  Fig.  1b.  The  exter¬ 
nal  airflow  conditions  and  the  liquid  rocket  nozzle 
exit  conditions  for  the  Titan  test  cases  are  pre¬ 
sented  in  Table  2.  The  Titan  three-dimensional  cal¬ 
culations  included  finite-rate  chemistry.  For  deter¬ 
mining  finite-rate  chemistry  effects,  the  chemical 
reaction  model  consisted  of  10  species  and  11 
reactions  to  represent  reactions  for  a  carbon, 
hydrogen,  oxygen  and  nitrogen-based  propellant 
system.  The  k-e  turbulenbe  model  was  used  for 
viscous  stress  approximation.  The  computation 
was  performed  on  a  single  node  of  a  four-proces¬ 
sor  SGI  Power  Challenge  R8000.  The  calculation 
executed  for  approximately  11,000  iterations,  and 
the  most  stressing  case  required  approximately 
four  weeks  of  CPU  time  to  converge  and  complete. 

Six,  three-dimensional  numerical  calculations  of 
the  twin-nozzle  configuration  were  obtained  at 
Mach  5.7  for  different  combinations  of  flow 
assumptions  and  approximations  in  an  attempt  to 


Table  1.  Parametric  Flow-field  Conditions 


Case 

Viscosity 

Inlet  Profile 

Chemistry 

Nozzle  Spacing 

1 

Laminar 

1-D 

Constant  y 

Narrow 

2 

Turbulent 

1-D 

Constant  y 

Narrow 

3 

Turbulent 

1-D 

Constant  y 

Wide 

4 

Turbulent 

1-D 

Finite  Rate 

Wide 

5 

Turbulent 

1-D 

Frozen 

Wide 

6 

Turbulent 

2-D 

Finite  Rate 

Wide 

3 
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a.  Schematic  of  the  3-D  plume 

Fig.  1 .  Geometry  representation  of  the  computational  domain. 


isolate  individual  influences  and  effects.  Cases  1  a  turbulent,  constant  gamma  approximation  with 

and  2  incorporated  laminar  and  turbulent  viscous  an  intranozzle  geometric  spacing  that  differed  from 

stress  models,  respectively,  and  were  computed  Cases  1  and  2.  Cases  4  and  5  were  obtained  with 

assuming  a  perfect  gas  equation  of  state.  Case  3  is  frozen  and  finite-rate  chemistry,  respectively. 
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Cases  1-5  all  assumed  uniform  nozzle  exit  proper¬ 
ties  as  the  GIFS  start  line  conditions.  The  final 
case  simulated  a  chemically  reacting  plume 
exhaust  using  radially  varying  profiles  for  the  initial 
conditions  at  the  nozzle  exit  plane  (calculated  via 
TDK)  to  define  the  starting  boundary  condition  for 
the  GIFS  plume  calculation. 

The  following  sections  will  discuss  the  individ¬ 
ual  influences  of  the  various  assumptions  for  phe¬ 
nomena  simulated  in  the  GIFS  model. 

Three-Dimensional  Effect 

In  an  earlier  study  by  one  of 
the  authors2  and  others,10  it  was 
shown  that  three-dimensional 
effects  are  important  and  should 
not  be  ignored  or  oversimplified  in 
modeling  efforts.  In  Ref.  2,  a  com¬ 
parison  of  the  two-dimensional 
axisymmetric  solution  with  the 
three-dimensional  twin-nozzle 
solution  was  accomplished  to 
determine  the  impact  of  the  single 
equivalent  nozzle  assumption. 

Figure  2  is  a  plot  of  static  pres¬ 
sure  contours  for  the  three- 
dimensional  solution.  The  plume 
expansion  shock,  barrel  shock, 
and  shock  reflection  are  clearly 
visible  at  1 16  meters  downstream 
of  the  nozzle  exit  plane  location. 

The  static  temperature  is 
increased  at  the  reflection  point  to  approximately 
half  the  value  of  the  total  temperature  of  the  flow.  A 
single  equivalent  nozzle  assumption  provides  valu¬ 
able  insight  concerning  overall  gross  qualitative 
assessments;  however,  for  detailed  studies  requir¬ 
ing  accurate  resolution  of  the  spatial  character  of 
multi  nozzle  plume  flows,  a  three-dimensional  cal¬ 
culation  is  required. 

Turbulence  Effect 

In  order  to  assess  the  influence  of  turbulence  in 
the  plume  flow-field  solution  resulting  from  twin 
nozzles,  two,  three-dimensional  calculations  were 
performed:  one  assuming  turbulent  flow  and  the 
other  assuming  laminar  flow.  As  would  be 


expected,  comparisons  of  the  flow-field  results  for 
the  turbulent  and  laminar  cases  indicates  decreas¬ 
ing  trends  in  plume  impingement  intensity,  plume 
mixing  rates,  and  a  smaller  plume  cross  section  for 
the  laminar  calculation. 

The  comparison  of  static  pressure  contours  in 
the  vertical,  horizontal  and  cross-planes  is  shown 
in  Fig.  3.  The  turbulent  and  the  laminar  solutions 
are  displayed  separated  centerline  and  indicated. 
Static  temperature  centerline  axial  profiles  con¬ 
trasting  the  turbulent  and  laminar  solutions  are 

Table  2.  Inflow  Conditions 


Fig.  2.  Pressure  contours  (3-D  calculation). 


Ambient  Conditions  at  47.6  km 

-  Tint  =  269  K  V  =  1 877.6  m/sec 

-  Cp/Cv  =  1.4  p  =  1 .388  x  1 0'3  kg/m3 

-P  =  108Pa  Mach  =  5.7 

-  Species  Concentrations 
(Mass  Fraction) 

-  N2  =  0.77  02  =  0.23 

Jet  Conditions 

-  T  =  1 920  K  V  =  2776.6  m/s 

-  P  =  92800  Pa  p  =  1 .68  x  1 0’1  kg/m3 

-  Mach  =  3.0 

-  Species  Concentrations 
(Mass  Fraction) 

CO  =  0.039  C02  =  0.1811  H2O=0.3496 

N2  =  0.414  NO  =  0.0109  OH  =  2.139e-3 

H2  =  3.13e-3  H=1.24e-4  O2  =  0.0  0  =  0.0 
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Z,  M  Z,  M 

Fig.  3.  Laminar  and  turbulent  pressure  comparisons. 


shown  in  Fig.  4.  These  results  indicate  that  the  bar¬ 
rel  shock  reflection  point  is  located  approximately  8 
meters  farther  downstream  for  the  turbulent  solu¬ 
tion.  Therefore,  the  spatial  characteristics  of  the 
plume  flow  field  and,  hence,  the  location  of  the 
radiation  centroid  can  be  influenced  by  the  viscous 
stress  model  included  in  the  solution.  In  the 
present  case  at  47.6  km  altitude,  the  effect  of  tur¬ 
bulence  is  not  particularly  strong.  However,  at 


Fig.  4.  Static  temperature  comparison  along 
centerline. 


lower  altitudes  and,  hence,  higher  Reynolds  num¬ 
bers,  turbulent  effects  will  be  more  significant. 

Intranozzle  Geometric  Spacing  Effect 

Changing  the  spacing  between  nozzles  should 
change  the  location  and  strength  of  the  plume/ 
plume  interactions.  In  order  to  explore  this  effect, 
perfect  gas,  three-dimensional  calculations  were 
obtained  at  two  different  intranozzle  geometric 
spacings  (narrow  versus  wide  spacing).  The  varia¬ 
tion  in  the  distance  between  nozzle  centerline  loca¬ 
tions  for  the  two  cases  was  approximately  1 5  per¬ 
cent,  with  the  first  case  assuming  the  widest  sepa¬ 
ration  distance.  Figure  5  compares  static  tempera¬ 
ture  contours  from  the  two  solutions.  These  results 
indicate  that  intranozzle  spacing  has  a  noticeable 
influence  on  the  flow-field  structure.  As  seen  in  Fig. 
5,  the  initial  plume  expansion  angle  is  larger  for  the 
wider  spacing  case  and  the  shock  reflection  loca¬ 
tion  is  farther  downstream.  A  centerline  axial  profile 
of  static  temperature  is  displayed  in  Fig.  6,  con¬ 
trasting  the  two  nozzle  spacing  solutions.  The  dif¬ 
ference  in  the  location  of  the  shock  reflect  point  is 
evident.  Figure  6  also  confirms  that  the  increased 
static  temperature  of  the  wide  spacing  case 
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Fig.  5.  Temperature  contour  of  perfect  gas  compari¬ 
son  nozzle  spacing  effect. 


extends  throughout  the  calculation  domain.  It  is 
concluded  that  the  intranozzle  spacing  affects  the 
plume  impingement  shock  location,  the  inviscid 
shock  structure,  and  the  shear  layer.  Thus,  nozzle 
spacing  can  significantly  impact  the  flow  field  and 
is  likely  to  influence  the  plume  emission. 

Kinetic  Chemistry  Effect 

A  frozen  solution  was  contrasted  with  a  finite- 
rate  reacting  solution  to  assess  the  significance  of 
chemistry  in  complex  plume  flow  fields.  A  compari¬ 
son  of  the  two  calculations  (Figs.  7-8,  respectively) 
indicates  that  there  is  an  effect  on  overall  plume 
structure  due  to  the  chemis¬ 
try  model.  The  increase  in 
static  temperature  attributed 
to  the  reflected  shock  is 
notably  different  between 
the  solutions.  The  barrel 
shock  reflection  for  the 
reacting  flow  solution  occurs 
farther  upstream  than  the 
frozen  case  and  results  in  a 
strong  normal  shock  Mach 
reflection/disc  behavior.  Fur¬ 
ther,  chemical  reactions 
occurring  in  the  plume 
impingement  and  shock 
reflection  region  result  in  dis¬ 
sociation  of  the  chemical 
species  and  also  recombi¬ 
nation  combustion  due  to 
mixing  with  the  air  in  the 


Fig.  6.  X-Y  plane  comparing  nozzle  spacing 
effect. 

region  between  the  nozzles.  This  phenomenon, 
which  was  not  observed  even  in  a  finite-rate  react¬ 
ing  single  equivalent  nozzle  solution,  is  due  to  mul¬ 
tiple  plume  interactions.  Figure  7  compares  temper¬ 
ature  contours  for  the  two  cases;  the  frozen  solution 
is  displayed  above  the  contour  centerline  and  the 
kinetic  solution  below.  The  presence  of  compres¬ 
sion-induced  combustion  in  the  downstream  shock 
region  is  evident  in  the  kinetic  case.  Figure  8  is  a 
centerline  axial  profile  of  static  temperature  extend¬ 
ing  from  the  nozzle  exit  plane  to  150  meters  down¬ 
stream.  This  comparison  indicates  that  the  chemis- 


100 

Axial  Distance,  M 

Fig.  7.  Static  temperature  contours  of  reacting  versus  frozen  flow. 
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try  model  has  an  effect  on  the  location  of  the  barrel 
shock  reflection  point.  The  reacting  chemistry  solu¬ 
tion  results  in  the  shock  reflection  point  occurring 
further  upstream  of  the  nozzle  exit  plane,  and  con¬ 
sequently,  an  overall  stronger  shock.  The  reacting 
chemistry  centerline  temperature  reaches  a  maxi¬ 
mum  value  of  over  2,000  K  behind  the  Mach  disc 
(Fig.  8).  For  the  frozen  case,  the  corresponding 
maximum  temperature  is  only  750  K. 

Initial  Boundary  Conditon  Profile  Effect 

In  a  previous  study  by  one  of  the  authors5  to 
assess  phenomena  affecting  scramjet  nozzle  per¬ 
formance,  it  was  determined  that  performance  is 
sensitive  to  variation  in  the  radial  profiles  assumed 
for  the  initial  boundary  condition  start  line  proper¬ 
ties.  Thus,  the  influence  of  the  nozzle  exit  plane  pro¬ 
file  variations  on  the  plume  flow  field  was  investi¬ 
gated.  In  one  case,  radial  exit  plane  profiles  were 
determined  via  the  TDK  model.  In  the  second 
case,  uniform  exit  plane  profiles  were  assumed. 
Plume  expansions  were  then  calculated  for  each 
case  assuming  fully  turbulent,  chemically  react¬ 
ing  flow.  A  comparison  of  static  temperature  con¬ 
tour  plot  and  cross  sections  at  axial  positions  of 
1  meter  and  7  meters  downstream  of  the  nozzle 
exit  are  shown  in  Fig.  9.  The  solution  displayed 
above  the  centerline  assumes  a  nonuniform  start 
line  profile.  The  solution  below  the  centerline 
assumes  a  one-dimensional  uniform  start  line 
condition.  Radial  profiles  of  static  temperature 
are  shown  at  axial  positions  equal  to  3,  10,  and 
20  meters  downstream  of  the  nozzle  exit  in 
Fig. 10.  A  comparison  of  the  spatial  plume  size 
indicates  that 
the  uniform 
starting  line 
results  in  a 
slightly  larger 
plume  expan¬ 
sion  region  due 
to  the  increased 
pressure  gradi¬ 
ent  between  the 
nozzle  exit  and 
the  free  stream. 

As  expected, 
differences 
between  the 


two  start  line  approximations  become  less  signifi¬ 
cant  as  the  flow  progresses  axially  downstream. 
Although  these  two  solutions  do  not  show  a  signifi¬ 
cant  effect  of  exit  profile  shape  on  the  plume  tem- 


Fig.  8.  Static  temperature  comparison  of  reacting 
versus  frozen. 


Fig.  9.  Temperature  contour  comparison  of  2-D  versus 
1-D  start  line. 


Fig.  10.  X-y  temperature  line  at  X  =  3, 10,  and  20  m. 
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perature,  it  is  expected  that  in  other  cases,  espe¬ 
cially  those  with  strong  afterburning,  that  larger 
effects  would  be  noted.  It  should  also  be  pointed  out 
that  the  aforementioned  solutions  were  obtained 
assuming  constant  oxidizer/fuel  (O/F)  ratio  across 
the  inflow  plane.  When  real  engine  effects  such  as 
fuel-film  cooling  and  injector  imperfections  are 
accounted  for,  O/F  can  vary  widely  across  the  noz¬ 
zle  exit  plane.  In  order  to  properly  account  for  the 
effect  of  such  mixture  ratio  variations  on  the  plume 
flow  field  and  resulting  base  heat  transfer  effects, 
radially  varying  exit  plane  profiles  must  be 
employed. 

Conclusions 

The  GIFS  code  has  been  significantly  modified 
to  improve  the  code  execution,  reliability,  and  user 
interface  for  three-dimensional  flow-field  calcula¬ 
tions.  A  modified  three-dimensional  GIFS  model 
was  applied  extensively  for  this  study  which 
included  six  different  three-dimensional  solutions. 
The  effects  of  nozzle  exit  nonuniformity,  intranozzle 
spacing,  turbulence,  and  finite-rate  chemistry  were 
evaluated  for  a  twin-nozzle/plume  propulsion  con¬ 
figuration.  Analysis  of  the  results  considering  these 
various  effects  leads  to  the  following  conclusions: 

1.  Nozzle  exit  start  line  assumptions,  nonuni¬ 
form  versus  uniform,  had  a  significant  effect  on  the 
immediate  plume  near  field.  The  difference 
between  the  two  start  line  assumptions  became 
less  significant  as  the  axial  distance  downstream 
of  the  nozzle  exit  plane  increased.  This  signifi¬ 
cance  of  the  observed  sensitivity  needs  additional 
investigation  to  elucidate  the  effect  of  the  starting 
conditions  on  the  plume  flow  field,  especially  in  the 
presence  of  strong  afterburning  and  missile  body/ 
base  flow-field  interactions. 

2.  Comparison  of  solution  results  contrasting 
laminar  and  turbulence  stress  models  indicates 
that  plume  flow-field  simulations  are  significantly 
influenced  by  the  turbulence  model.  A  laminar 
approximation  is  not  appropriate  for  use  in  low-  to 
moderate-altitude  plume  flow-field  analysis. 

3.  The  chemistry  assumption  (frozen  versus 
reacting)  did  not  influence  the  overall  plume  struc¬ 


ture  in  this  case.  However,  the  chemically  react¬ 
ing  solution  produced  a  strong  Mach  reflection 
behavior  and  combustion  in  the  vicinity  of  the 
Mach  reflection  zone  which  was  not  predicted  by 
the  frozen  results. 

4.  The  intranozzle  spacing  distance  has  a  sig¬ 
nificant  impact  on  the  barrel  shock  reflection  loca¬ 
tion,  plume/plume  impingement  shock  location, 
and  shear  layer. 

5.  A  comparison  of  the  two-dimensional  and 
three-dimensional  cases  indicates  that  the  three- 
dimensional  effects  are  important  in  the  near-field 
plume  and  diminish  as  the  axial  distance  extends 
further  downstream  from  the  nozzle  exit  plane.  The 
single,  equivalent  nozzle  approach  should  not  be 
used  to  describe  plume  near-field  flow  characteris¬ 
tics  where  three-dimensional  effects  predominate, 
and  in  the  instances  where  3-D  features  are 
required  as  part  of  the  flow-field  description. 

6.  These  analyses  indicate  that  two-dimen¬ 
sional  nonequilibrium  analysis  tools  can  provide 
some  insight  concerning  overall  gross  qualitative 
assessments  of  multiple  plume  flow-field  phenom¬ 
ena.  However,  for  detailed  studies  of  complex 
flow-field  phenomena,  a  more  sophisticated  three- 
dimensional  calculation  is  required. 

Furthermore,  this  study  is  intended  to  assess, 
understand,  and  quantify  phenomena  of  the  plume 
physics  and  to  identify  where  simplified  models  are 
appropriate  for  use  without  introducing  any  signifi¬ 
cant  error  to  the  simulated  flow  field  and  conclu¬ 
sions  which  might  be  deduced  from  analysis  of  the 
flow-field  simulation. 

Future  extensions  to  the  study  presented  in  this 
paper  will  explore  three-dimensional  vehicle/base 
interactions  with  the  plume  flow.  Also,  to  promote 
optimal  utilization  of  computer  resources,  hybrid 
Navier-Stokes/Parabolized  Navier-Stokes  methods 
are  being  investigated. 
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